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Abstract 

While Hermiticity of a time-independent Hamiltonian leads to unitary time evolution, in and of 
itself, the requirement of Hermiticity is only sufficient for unitary time evolution. Here we provide 
conditions that are both necessary and sufficient. Specifically, we show that VT symmetry of a 
time-independent Hamiltonian, or equivalently, reality of the secular equation that determines its 
eigenvalues, is both necessary and sufficient for unitary time evolution. For any 'PT-symmetric 
Hamiltonian H there always exists an operator V that relates H to its Hermitian adjoint according 
to VHV^ 1 = W. Hilbert space norms (i/ji\V\ip2) constructed with this V are always preserved 
in time. With the energy eigenvalues of a real secular equation either being real or appearing in 
complex conjugate pairs, we thus establish the unitarity of time evolution for both stationary and 
non-stationary systems. In addition, we show that such systems obey causality. 
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I. UNITARY TIME EVOLUTION FOR NON-HERMITIAN HAMILTONIANS 



If for a time-independent Hamiltonian H one constructs a time evolution operator of the 
form U(t) = e~ tHt , then if H is Hermitian, the adjoint operator W(t) is given by lP(t) = e lHt 
and the unitarity of U(t) immediately follows. Moreover, if one introduces right-eigenvectors 
of H that obey 

il\Ri(t)) = H\Ri(t)), (1) 

then for Hermitian Hamiltonians the Dirac norm (Ri(t)\Ri(t)) is immediately preserved in 
time. However, if the Hamiltonian is not Hermitian, conjugation of (CQ) leads to 

- i^imi = im\H\ (2) 

with U\t) now being given by U\t) = e lHtt , and with the Dirac norm then not being time 
independent. Despite this, one cannot conclude that the time evolution associated with a 
non- Hermitian Hamiltonian is not unitary. Rather, all one can conclude is that the standard 
derivation of unitarity for Hermitian Hamiltonians and the appropriateness of the Dirac norm 
do not apply when the Hamiltonian is not Hermitian. While the above analysis shows that 
Hermiticity implies unitarity, one cannot conclude that lack of Hermiticity implies lack of 
unitarity. And indeed, a whole class of non-Hermitian Hamiltonians has been identified 
in which time evolution is unitary, namely those in which the Hamiltonian possesses a 
VT symmetry and has a set of energy eigenvalues that are all real (see the review of [l| 
and references therein). The purpose of this paper is to identify all possible classes of 
Hamiltonians for which time evolution is unitary, with Hermitian Hamiltonians being just 
one of the possible classes. 

The key to finding the general rule is to recognize that if \Ri(t)) is an eigenstate of 
H with some general complex energy E\ = Ef + iEf, then on introducing some general 
time-independent operator V, the ^/-dependent inner product (Rj(t)\V\Ri(t)) will be given 
by 

(Rjit^VlRiit)) = e-^-^+^'+^iRjit = 0)\V\Ri{t = 0)). (3) 

This inner product will be time independent for energies that obey the two conditions 

Ef - Ef = 0, E{ + E] = 0, (4) 

and will be time independent for states with energies that do not obey these conditions 
if an operator V can be found so that for all of those states the matrix element (Rj(t = 



0)|V|i2i(i = 0)) vanishes. The eigenstates of Hermitian Hamiltonians immediately obey 
these conditions with E\ = 0, Ej = and V = I. However in the general case one needs a 
non-trivial V, and a thus non-Dirac norm. 

To find the general conditions under which the needed V will in fact exist, we note that 
from (JTJ) and we can evaluate the time derivative of (Rj(t)\V\Ri(t)) to obtain 

i±(W)\v\W)) = (RM(vh - H*v)\RiV))- (5) 

Thus if, just as is needed to establish unitarity for Hermitian Hamiltonians, we require that 
the set of states \Ri(t)) be complete for the general H, we see that all matrix elements 
(Rj(t)\V\Ri(t)) will be time independent if V obeys the relation 

VHV' 1 = (6) 

as an operator identity. The needed V is thus the one that effects the transformation from 
H to its Hermitian adjoint H'. Conversely, if a Hamiltonian has the property that H and 
can be related by an operator V as in ([6]), it follows from (jHJ) that all matrix elements 
(Rj(t)\V\Ri(t)) constructed with this V are time independent. The condition that there 
exist an operator V that effects VHV^ 1 = H* is thus a necessary and sufficient condition 
for unitary time evolution. We note that as defined in ([6]) the operator V depends on the 
Hamiltonian. Thus unlike the standard Dirac norm (Rj(t)\Ri(t)), the (Rj(t)\V\Ri(t)) norm 
cannot be preassigned, with it being dynamically determined by the Hamiltonian itself. 
When there does exist an operator V for which ([6]) is satisfied, we can rewrite (j2J) as 

-i±(R t {t)\V=(R i (t)\VH. (7) 

In such a situation the state (Li(t)\ = {Ri(t)\V is a left-eigenvector of H as it obeys 

-i^(I H (t)\ = (L i (t)\H 1 (8) 

with unitary time evolution thus being associated with Hamiltonians that have complete 
sets of both left- and right-eigenvectors. In terms of the left-eigenvectors the (Rj(t)\V\Ri(t)) 
inner product can be written as (Lj(t)\Ri(t)), with the inner product thus being the overlap 
of the left- and right-eigenvectors of the Hamiltonian, an overlap that is manifestly time 
independent because of (CQ) and (jSJ). The general rule then is to define the inner product as 
the overlap of the left- and right-eigenvectors, and so we need to ascertain the conditions 
on H for which a V operator that satisfies ([6]) (and thus the left-eigenvectors that it then 
generates) will necessarily exist. 
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II. VT SYMMETRY AND THE VHV" 1 = CONDITION 

When H and W are related according to the similarity transform VHV^ 1 = H\ the 
determinant condition for the eigenvalues of H obeys det(if — XI) = detiV^H^V — XI) = 
det(iJ t - XI) = det((#*) TR - XI) = det{H* - XI) = 0, where TR denotes transpose. 
Consequently, H and H* possess the same set of eigenvalues, with the secular equation 
det(H — XI) = thus being real. With the energies of a real secular equation either being 
real or appearing in complex conjugate pairs, one can immediately satisfy (jlj) non-trivially, 
just as required. Conversely, if the secular equation that determines the energy eigenvalues 
is real, then H and will possess the same set of eigenvalues and thus be related by an 
isospectral similarity transform. Reality of the secular equation is thus both necessary and 
sufficient for the existence of a V operator that obeys §6§. Reality of the secular equation is 
thus both necessary and sufficient for unitary time evolution. 

The reality requirement on a secular equation is a condition that is encountered in VT- 
symmetric theories (V and T denote parity and time reversal). Specifically, it was noted in 
2| that if a Hamiltonian is VT symmetric, the secular equation would be real. Moreover, it 
is not specifically VT invariance itself that is required for the proof. All one needs is that H 
be invariant under the action of the product of a discrete linear operator and a discrete anti- 
linear operator. For a time-independent H, time reversal reduces to complex conjugation 
K times a matrix that acts in the space of the eigenvectors of H, with VT then acting as 
AK where A is a matrix. The essence of the proof of the reality of the secular equation is 
to note that if a Hamiltonian commutes with AK, the determinant condition then obeys 
det(# - XI) = detiAKHKA- 1 - XI) = det(KHK - XI) = det(H* - XI) = 0, with reality 
of the secular equation immediately following. Discussion of the reality of secular equations 
dates back to Wigner {3J], who considered invariance of a Hamiltonian under time reversal 
alone. Modern interest in the reality of secular equations was triggered by the surprising 
discovery that the non-Hermitian Hamiltonian H = p^/2m + ix 3 had a completely real 
eigenspectrum 0, S], an outcome that was traced in [5J to the fact that while not time 
reversal invariant, H was nonetheless invariant under the VT product as T transforms i 
into —i while V transforms x into —x. The utility of VT invariance stems from the fact that 
complex potentials that are not T invariant might still be invariant with respect to a VT 
product where V is an appropriate linear operator, and thus yield a real secular equation. 
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In |2| it was shown that if a Hamiltonian is VT symmetric then its secular equation is 
real. More recently, the converse has been proven [(J, namely that if a Hamiltonian has a 
secular equation that is real, then one can always find appropriate V and T operators under 
which the Hamiltonian is VT symmetric. Consequently, VT symmetry of a Hamiltonian is 
both necessary and sufficient for reality of its secular equation, and thus VT symmetry of a 
Hamiltonian is both necessary and sufficient for unitary time evolution. The power of this 
result is that one can determine whether or not a Hamiltonian can generate unitary time 
evolution entirely by examining its symmetry structure under the action of the product of 
a discrete linear operator and a discrete anti-linear operator. 

Within the family of Hamiltonians that have real secular equations there are three cases 
of interest: (i) the energies are all real and the energy eigenvectors are complete, (ii) the 
energy eigenvectors are complete but the energies include one or more complex conjugate 
pairs, and (iii) the energies are real but the energy eigenvectors do not form a complete set 
and the Hamiltonian is in non-diagonalizable, Jordan block form. A possible fourth case, 
namely that the energy eigenvectors are not complete and the energies include one or more 
complex conjugate pairs, is not allowed since the energies of a Jordan-block matrix are all 
equal (they all share a single common eigenvector) and could thus not contain any complex 
conjugate pairs |6|. Moreover, in 6| criteria were also provided for determining to which of 
these three cases of interest any given PT-symmetric Hamiltonian belongs. Specifically, in 
the general analysis of PT-symmetric Hamiltonians it has been found to be very convenient 
[l| to introduce a linear operator, the C operator, that obeys the two conditions: [C, H] = 0, 
C 2 = 1. In [6] it was shown that if every C that non-trivially obeys these two conditions 
commutes with VT then we are in case (i), that if there exists at least one C that does 
not commute with VT then we are in case (ii), and if there is no C that obeys these two 
conditions at all then we are in case (iii). Thus as well as ascertain the unitarity structure 
of a Hamiltonian via a symmetry condition (i.e. whether or not H commutes with VT), we 
can even determine which realization of the energy spectrum we are in by examining the 
structure of a commutator (viz. that of C with VT) . We now examine how the unitarity of 
time evolution works in each of these three cases in detail. 
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III. UNITARITY WHEN ENERGIES ARE ALL REAL AND COMPLETE 



For the case of a PT-invariant Hamiltonian whose energy spectrum is both real and 
complete, it can be shown that the operator V can always be written as e ^ 

where Q is a Hermitian operator, and that the operator e~® itself can always be written 
as e"^ = PC where C is the C operator discussed above. The utility of the Q operator is 
that its Hermiticity ensures that the operator H = e~®/ 2 He®/ 2 is Hermitian [lllsj]- If for 
H we introduce a standard set of right-eigenstates \n l ), the right-eigenstates of H are thus 
related to those of H by the mapping |rtj) = e^^ 2 \hi). Since Q is Hermitian and e~®/ 2 is a 
similarity transform but not a unitary one, the states (rii\ = (n^e^l 2 are not left-eigenstates 
of H. Rather, it is the states {n^e"^ that are the left-eigenstates. Moreover, since the 
eigenstates of the Hermitian H obey the standard time-independent Dirac orthonormality 
condition (ji l \n^) = S^j, then by the mapping the eigenstates of H obey the time-independent 
condition {n l \e~^\n^) = d^j. With V = e~®, we recognize the inner product {n l \e~^\n^) as 
being none other than the time- independent ^-operator norm introduced above [9j. 

The essence of the above analysis is the realization that when a Hamiltonian H has an 
energy spectrum that is both real and complete, the Hamiltonian is either already Hermitian 
or can be brought to a Hermitian form H by a similarity transformation. However, if a 
Hamiltonian H is not Hermitian the needed similarity transformation will not be unitary and 
the basis states of H will not be orthonormal with respect to the Dirac inner product. Rather, 
they will form a skew basis, with it being the similarity transform that puts the basis into 
an orthonormal form. PT-symmetric Hamiltonians that are not Hermitian but nonetheless 
have a real and complete eigenspectrum are thus Hermitian in disguise. Consequently, they 
are unitary in disguise too, with it being the ^-operator norm that is the appropriate norm 
for the Hamiltonian. 

It is instructive to recast the above discussion in the language of the S'-matrix. For 
the above \hi) basis first, we can define complete sets of in and out states that obey the 
standard closure relations J2i I^LX^iJ = I an d J2i \n l ut) i^outl = I- Given these relations 
one can show that the scattering operator defined by S = J2i \^l n ) (^out\ an d thus by = 
J2i \nlut) falnl immediately obeys the unitarity relation SS^ = I. 

Since the in and out states of H can be mapped into the in and out states of H by the 
mapping, an S'-matrix operator in the |rij) basis can thus be defined via the mapping as 
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S = e~®/ 2 Se®/ 2 , with the adjoint being given via = e®/ 2 S^e~"®/ 2 since the operator Q 
is Hermitian. As constructed, the operator S can be written as S = e^ 2 J2i I^L) (^otJ e ~^ 2 \ 
and from the mapping thus takes the form S = J2i \ n in) i n out\ e ~ Q i n the \rii) basis. In 
addition, the SS^ = I relation obliges S to obey e -Q/ 2 Se Q l 2 e Q l 2 S^e~ Q l 2 = I, i.e. to obey 

Se Q S*e~ Q = I, e- Q Se Q S ] = I. (9) 

It is thus 02]) that is the form that unitarity of a scattering process takes in a skew basis. 

While we have derived (jHJ) by transforming to the Hermitian basis and then transforming 
back, it is possible to derive (Q directly in the \nj) basis without needing to make any 
reference to the tilde basis at all. Moreover, there is a more general derivation that does not 
even require us to restrict to Hermitian V. Specifically, we note that for energy eigenstates, 
be they real or complex, the time independence of (Rj(t)\V\Ri(t)) as required by ([51) and 
(jSJ) entails that states whose energies do not obey (BJ have to obey (Rj(t = 0)\V\Ri(t = 
0)) = 0. Consequently, such states are orthogonal with respect to the V^-operator norm, 
and on normalizing the states appropriately, for all energy eigenstates we can thus set 
(Rj(t)\V\Ri(t)) = Sij. With this orthonormality relation we obtain orthonormality and 
closure relations for the in and out states of the form 

(K n \V\ni n ) = (rijVVL} = Sij, K ut \V\ni ut ) = (nLl^VD = ^ 

EW(^ = E yt W(<l = ^ T,\<ut)(<ut\v = T,v j Kut)K ut \ = i-m 

i i i i 

If we now define an operator S = J2% \ n l n ) i n out\V > from (fTD]) we obtain 

SV-WV = £ InDKJVV-^V^nUinllV = £ K>«JV = /. (11) 

i j i 

Consequently, for PT-invariant Hamiltonians unitarity of the scattering process thus takes 
the equivalent forms 

SV- l S ] V = I, VSV^S* = I. (12) 

As a simple check on (fT2|) . we note that if S is written as the evolution operator S = e~ lHt , 
then because = VHV -1 , it follows that VSV' 1 is given as e~ tHlt . Then with Sft = e +iiTt *, 
(TT2"j) immediately follows. The relation given in (I12j) is thus the unitarity relation we want. 
While we shall discuss the complex energy case in more detail below, as we had already 
noted, our derivation of (I12j) holds even if some energies are complex. The issue here is 
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not whether or not we can prepare stable in and out states with them, but whether or not 
we need them for the closure relations given in ([TO]) . In other words, suppose we have a 
PT-invariant Hamiltonian with N real energy eigenvalues and M complex pairs. For such 
a situation, unitarity in the form given in f[T2"j) will only follow if in the closure relations in 
(HOP we sum over iV + 2M states (i.e. over iV + 2M channels), and not over N states alone. 

As well as provide a formal derivation of our results, it is instructive to see how our ideas 
work in a simple model. Since we can always diagonalize any Hamiltonian that possesses a 
complete set of eigenstates, and since complex congugate pairs would form a two-dimensional 
subblock in the diagonal basis, it will suffice to discuss a two-dimensional model. We thus 
follow [l| and introduce the Hamiltonian (the parameters r, s and 9 are real) 

/ r cos 9 + ir sin 9 s \ 

H = r cos 9ao + ir sin 0cr 3 + so\ = . (13) 

\ s r cos 9 — ir sin 9 J 

While not Hermitian, this H is VT invariant under V = a\, T = K and thus has a real 
secular equation with energies being given by E± = r cos 6 ± (s 2 — r 2 sin 2 6) 1 / 2 . In the region 
where s 2 — r 2 sin 2 9 is positive the two energy eigenvalues are real and distinct. In this region 
we can define a Hermitian operator e~®l 2 according to 

to/2 /l + sina\ 1/2 /l-sina\ 1/2 TO 1 , cosa 



2 sin a J V 2 sin a J sin a sin a 

where sin a = +(s 2 — r 2 sin 2 9) 1 ! 2 / s, cos a = r sin 9/s. Using the e~®/ 2 operator we construct 

H = e' Q ' 2 He Q ' 2 = r cos 6a + a^s 2 - r 2 sin 2 9) 1 ' 2 , (15) 

and confirm that H is Hermitian when s 2 — r 2 sin 2 9 is positive. Similarly we construct 
e~®He® and confirm that it is equal to W. Finally, we construct an operator C = e®P, 
viz. 

1 

C— (di + i cosaa"3) , (16) 

sin a 

and confirm that it not only obeys [C, H] = and C 2 = I, it also obeys [C, VT] = 0, just as 
it must since the energies are real. 

In the region where s 2 — r 2 sin 2 9 is positive, eigenvector solutions that obey idu±/dt = 
Hu± = E±u± are given as 

g— i(rcos0+^)tgi-7r/4 / g — ' a / 2 \ ^--i(r cos 9 — [i)t ^iir / A / ^g* Q / 2 \ 

U+ = (2sina) 1 /2 I _ K »/2 J ' M - = (2sina) 1 /2 1 e ~ia/ 2 J ' ( 17 ) 



S 



where /i = + (s 2 — r 2 sin 2 #) 1//2 . For the states in (|17[) the e ^ norms obey the time- 
independent and manifestly unitary orthonormality and closure relations 

u t ±e~ Q u± = +1, v)±e- Q u T = 0, u + u\e~ Q + u^ule~ Q = J, (18) 

just as required. It is thus the V = e~ Q norm that it is the appropriate one for the problem. 

IV. UNITARITY WHEN ENERGIES ARE NOT ALL REAL BUT COMPLETE 

To determine what happens when energies are complex, it is instructive to study the 
two-dimensional model in the regime where s 2 — r 2 sin 2 9 is negative. Now the energies are 
given by the complex conjugate pair E± = rcosO ± iv where v = +(r 2 sin 2 # — s 2 ) 1 / 2 . In 
this region we can define a now non-Hermitian operator e~^/ 2 according to 

t q/2 ( 1 + * sinh /A ( 1 — i sinh (3 \ 1 ^ 2 ^.q 1 cosh/3 

\ 2i sinh (3 J 2 \ 2i sinh (3 J ' i sinh (3 2 i sinh (3 ' 

(19) 

where sinh/? = +(r 2 sin 2 # — s 2 ) 1 ^ 2 /s, cosh/? = rsin#/s. Using this e~®/ 2 operator we 
construct 

H = e- Q/2 He Q/2 = r cos 6a + ia x (r 2 sin 2 9 - s 2 ) 1/2 , (20) 

and see that H is not Hermitian when s 2 — r 2 sin 2 9 is negative. Despite this, we construct 
e~"®He® with this non-Hermitian Q and confirm that it nonetheless is equal to H\ just as 
it must be. Finally, we construct an operator C = e®P, viz. 

C=^—-(a 1 + icoshPa 3 ) (21) 
i sinh [3 

and confirm that it while it obeys [C, H] = and C 2 = I, because of its overall factor of i it 
does not obey [C,VT] = 0, as must be the case since the energies are not real. 

In the region where s 2 — r 2 sin 2 # is negative, eigenvector solutions that obey idu±/dt = 
Hu± = E±u± can be constructed by setting /i = iv, a = i(3 in f|T7|) . and are given as 

g— ir cos 0t+ut I C^/ 2 \ ^—irca&Qt—vt /ic~@l 2 \ 

U+ = (2sinh/?) 1 /2 [ _ ie -P/2 ) ' U ~ = (2sinh/?) 1 /2 [ e p/2 ) ' (22) 

For the states in (1221) the e~® norms obey the time-independent orthogonality and closure 
relations 

v}±e~ Q u± = 0, u^e~~ Q u + = +1, u\_e~ Q u_ = -1, 

u + ule- Q - u^u{e- Q = I. (23) 



With the adjoints of the states in (1221) having time dependences of the form u\ ~ e vr cos et+vt ^ 

that in (123]) those overlaps that are zero have to vanish through the 
matrix structure of the overlaps since the time dependences in u\-e~®u± do not cancel 
each other once v is non-zero. For the non-zero overlaps the time dependences in u^e~^u± 
do precisely cancel each other, just as required by PJ. With the u + and w_ states in (122]) 
obeying VTu + = u_, VTu_ = u + (T effects t — ► — £ as well as z — * — i), the non- vanishing 
overlaps in (1231) represent transitions between energy eigenstates and their VT partners jlO ]. 

The most notable aspect of the overlaps is that the u+e~®u_ overlap is negative. How- 
ever, this does not herald the presence of a negative metric since the overlap is a transition 
matrix element between different states and not the overlap of a state with its own ad- 
joint. Moreover, rather than being an indicator of a possible loss of unitarity, this negative 
sign is actually needed for unitarity. Specifically, if we construct a propagator D(E) for 
above Hamiltonian, it will have the form of matrix elements of the states divided by energy 
denominators: 

D{E) = U \ 6 ' QU+ + U 'f Qu - , (24) 
1 ' E-(E R + iEj) E - (E R — iEj) 1 ; 

and thus evaluate to 

D{E) = E-(E R + ify) " E — (E R — iEj) = (E - E R ) 2 + E]' (25) 
Other than the overall factor of 2, the imaginary part of this propagator is the same as that 
of a standard Breit-Wigner 

n fjp\ 1 E-E R + iEj 

Dbw(E) = E-iExTW) = (E — E R ) 2 + E] ' (26) 

The negative sign of u\_e~®U- thus compensates for the negative sign of the imaginary part 
of the energy of the w_ state, with its contribution to the imaginary part of D(E) then being 



positive rather than negative 



While it is nice to see that D(E) recovers the standard Breit-Wigner form, our results 
actually go further than that as they provide some justification for the use of the standard 
D^{E) in the first place. Specifically, if one works with a general complex potential, i.e. 
with one for which (j4j) is not obeyed, one would not get unitary time evolution. If for instance 
we work with a complex potential that has just one energy eigenstate with wave function 
ip = e -i' E R t + E i t ^ the probability would be given as ip*ip = e~ 2Ejt and not be preserved in time, 
and would even grow in time if Ej is negative. Hence, if we wish to describe absorption or 
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decay processes in a time-preserving manner, we must augment this state with its complex 
conjugate partner and include the contribution of the partner to D(E). Our approach is thus 
well-suited to a scattering process in which the in and out states long before and long after 
a collision are real eigenvalue eigenstates of a PT-invariant Hamiltonian H , with the states 
scattering through some short-range PT-invariant interaction Hint f° r which H + i?TNT 
has complex eigenvalues. In such a situation the interaction would generate two poles at 
Er ± iEj that would be on the same Riemann sheet and combine just as in (125]) 121 ] . 

The fact that the two complex conjugate poles are on the same Riemann sheet has an 
interesting consequence: it enables the two poles to be located within the same contour in the 
complex energy plane rather than be in different contours. In consequence, the theory will be 
causal. Specifically, with two complex conjugate poles having opposite imaginary parts, one 
of the poles would have to appear above the real energy axis and one below. Now when one 
only has to deal with poles with real energies (such as the D(E) = l/E + + 1/E_ propagator 
associated with ( fl8l) ). to construct a causal, retarded propagator, one locates all the poles 
slightly below the real axis, with a contour integration then giving no pole contribution if 
one closes the contour in the upper-half plane. If one simply repeats this prescription when 
there is a pole in the upper-half plane, on closing the contour above the real axis one then 
gets a non-zero pole contribution and causality violations result. However, as discussed 
for instance in [13j, for complex conjugate pairs of poles one can use an alternate contour 
instead. Specifically, one can deform the contour to pass above any pole in the upper-half 
plane. Then when one closes this deformed contour in the upper-half plane, one still gets no 
pole contribution and the propagator is causal, and when one closes the deformed contour in 
the lower-half plane, both complex poles contribute, just as in (1251) . For VT theories then, 
to maintain causality, the prescription is to make the choice of contour dynamics dependent 
(just as is done with the PT-theory Hilbert space metric), so that as interactions cause 
poles to move off the real axis the contour must move along with them so that all the poles 

n 

remain below it, a meaningful procedure if all the poles are on the same Riemann sheet |14| . 

V. UNITARITY WHEN ENERGIES ALL REAL BUT INCOMPLETE 

When a Hamiltonian is diagonalizable, one does not need to distinguish between the 
eigenvalue solutions to the secular equation det(H — XI) = and the eigenvalue solutions 



11 



to Hip = Eip since the two sets of solutions coincide. However, when a Hamiltonian is 
not diagonalizable, the number of eigensolutions to Hip = Eip is less than the number of 
eigenvalue solutions to det(H — XI) = 0, and the spectrum of H is incomplete. A typical 
example of a non-diagonalizable, Jordan-block, matrix is the two-dimensional matrix M 

M =C (27) 

Its secular equation \M — XI\ = has two solutions for A, both of which are equal to a 
(and incidentally both solutions are real when a is real even though an M with real a is not 
Hermitian), but M has only one right-eigenvector and equally only one left-eigenvector, viz. 

(o !)C)= a C)- 0(0 i)=(o D- (28) 

Despite the lack of completeness of its eigenstates, if the secular equation det(H-XI) = 
of a Hamiltonian is real, H and its adjoint W will still be related by ([6]), i.e. still be related 
VHV~ l = W for some operator V. For the example of (1281) the secular equation will 



be real when a is real, and V will then be given by V — G\ [15j]. More generally, for a 
non-diagonalizable matrix of arbitrary dimensionality, Jordan showed that via a sequence of 
similarity transformations any such matrix can always be brought to the Jordan canonical 
form in which all of the elements on the diagonal are equal to each other, and the only non- 
zero off-diagonal elements are all real and all lie on the superdiagonal immediately above the 
diagonal. With the non-zero elements on the superdiagonal not contributing to the secular 
equation since all the other non-diagonal elements of the matrix are zero, the elements on 
the diagonal of a Jordan block matrix are the eigenvalue solutions to the secular equation. 
Since all the eigenvalues of a Jordan-block matrix have to be equal to each other, the reality 
of the secular equation of a Jordan-block matrix requires that all of its solutions be real, 
since any complex conjugate pair of solutions could not be equal to each other. Thus for 
any Jordan-block matrix whose secular equation is real, we see that all the elements of the 
matrix are real, with the Hermitian adjoint of the matrix thus being equal to its transpose. 
Since a transposition matrix always exists in any dimension (an explicit construction is given 
in 6|), the transposition matrix serves as V. Thus for any Jordan-block matrix with a real 
secular equation, one can always find an operator V that effects VHV^ 1 = W . 

Moreover, in the Jordan block case the converse also holds. Specifically if one can find a 
V operator that effects VHV~ l = W , the secular equation will be real since the invariance 
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of a determinant under similarity transforms is insensitive to Jordan-block structures. Thus 
even in the Jordan block case VT invariance is necessary and sufficient for the existence of 
a V operator that effects VHV -1 = H'. 

With (jSJ) holding for any states that obey ([T]) and (J2J), flSJ) will hold for non- stationary as 
well as stationary solutions to (UJ. Now in the Jordan block case the lack of energy eigenstates 
is compensated for (see e.g. [16|) by the presence of non- stationary solutions that typically 
behave as powers of t. For every missing stationary solution one finds a power-behaved 
in t solution The set of all stationary plus non-stationary solutions combined is thus 
complete [li|. Since the stationary plus non-stationary solutions are complete, we can use 
(J5J) to infer that in the Jordan block case time evolution of wave packets will be unitary if 
and only if one can find a V operator that effects VHV^ 1 = W. Thus, as noted in |l6| . 
even though the non-stationary solutions are power-behaved in t, packets containing them 
and the stationary solutions will still evolve with a probability that is preserved in time fl9 ]. 
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